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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  new  and  emerging  analysis  and  forecast  capabilities 
via  improvements  to  environmental  inputs  to  electro-magnetic  (EM)  tactical  applications  through 
Unmanned  Aerial  Systems  (UAS)  enhanced  high-resolution  mesoscale  data  assimilation  and  modeling. 
These  capabilities  will  be  achieved  through  the  formulation,  testing,  evaluation,  and  demonstration  of 
an  environmental  data  assimilation  and  prediction  system  designed  to  quantify  the  UAS  observation 
parameters,  characteristics  (e.g.,  vertical  resolution,  rate,  etc.),  measurement  strategies  and 
assimilation  approaches  that  provide  the  greatest  benefit  for  high-resolution  (grid  increment  of  5-km  or 
less)  atmospheric  tactical-scale  environmental  forecasts.  The  goal  is  to  focus  on  EM  tactical 
applications  and  metrics  in  evaluating  the  performance  of  the  system.  This  development  will  leverage 
current  and  emerging  capabilities  for:  i)  environmental  data  assimilation,  ii)  adjoint-based  sensitivity 
and  observation  impact,  iii)  air-ocean  coupled  modeling,  and  iv)  EM  tactical  decision  aid  systems. 
Observational  results  from  field  programs  such  as  Trident  Warrior  2013  (TW13)  will  also  be  leveraged 
to  build,  test,  and  transition  the  prototype  end-to-end  (observation  to  tactical  decision  aid  (TDA)) 
demonstration  system.  All  meteorological  inputs  to  the  demonstration  system  will  come  from  UAS 
platforms  and  sensors  consistent  with  the  present  state-of-the-art;  however,  no  specific  operational 
UAS  system  will  be  targeted,  transitioned  or  modified  as  part  of  this  effort. 

OBJECTIVE 

The  objective  is  to  develop,  validate,  and  transition  new  METOC  capabilities  within  a  high-resolution 
modeling  system  using  advanced  data  assimilation  of  UAS  data  in  order  to  achieve  superior  short  and 
medium  term  forecasts  of  high-impact  weather  and  EM-related  parameters.  We  will  use  three  crucial 
model  components:  i)  the  coupled  atmospheric  and  oceanic  components  of  the  Coupled 
Ocean/Atmosphere  Mesoscale  Prediction  System  (COAMPS®1),  ii)  the  3D-variational  (3DVAR) 
analysis  schemes  NAVDAS  (NRL  Atmospheric  Variational  Data  Assimilation  System)  for  the 
atmosphere  and  NCODA  (Navy  Coupled  Ocean  Data  Assimilation)  for  the  ocean  and  wave  models, 
and  iii)  the  COAMPS  Adjoint  observation  impact  system.  COAMPS-OS®1  will  be  used  as  the  software 
interface  for  configuring  COAMPS  and  the  other  system  components  to  facilitate  high-resolution 
coupled  mesoscale  forecasts.  The  COAMPS  system  includes  an  ocean  circulation  model,  the  Navy 


1  CO  AMPS®  and  COAMPS-OS®  are  registered  trademarks  of  the  Naval  Research  Laboratory. 
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Coastal  Ocean  Model  (NCOM),  which  fully  interacts  with  the  nonhydrostatic  atmospheric  module.  In 
addition,  two  different  wave  models  are  currently  being  implemented  and  tested:  Simulating  Waves 
Nearshore  (SWAN)  and  Wavewatch  III  (WWIII)  which  will  also  be  available  for  testing.  NAVDAS  is 
the  Navy’s  state-of-the-art  analyses  system.  The  adjoint  of  the  CO  AMPS  forecast  model  and 
NAVDAS  data  assimilation  systems  have  already  been  developed  and  linked  together  to  provide  a 
quantitative  measure  of  the  impact  of  individual  observation  types  (e.g.,  radiosonde,  aircraft,  UAS, 
etc.)  on  the  COAMPS  forecast  error.  This  observation  impact  capability  has  only  recently  become 
available  for  COAMPS  and  is  one  of  the  first  such  observation  impact  capabilities  for  a  mesoscale 
modeling  system  developed  anywhere. 

Several  key  scientific  objectives  are  identified  in  relation  to  the  impact  of  the  environment  on  EM 
propagation  characteristics,  particularly  in  the  coastal  marine  atmospheric  boundary  layer  (MABL):  i) 
the  time-varying  3-D  structure  of  the  MABL  and  its  interaction  and  impact  on  EM  propagation 
characteristics  and  ducting  of  EM  energy,  ii)  the  response  of  the  MABL  to  the  structure  of  the 
underlying  SST,  including  SST  gradients,  and  the  impact  on  evaporative  duct  height  calculations,  on 
MABL  refractive  layers,  and  on  larger-scale  surface  and  MABL  ducting  phenomena,  and  iii)  the  most 
effective  UAS  observing  strategies  from  an  atmospheric  model  data  assimilation  perspective  in  order 
to  accurately  predict  mesoscale  weather  and  EM  propagation  conditions. 

APPROACH 

Our  approach  is  to  leverage  prior  investments  in  NRL-developed  NAVDAS  (including  aircraft  data 
assimilation),  COAMPS,  adjoint,  and  EM  TDA  technologies  to  produce  a  UAS-enabled,  integrated, 
and  automated  on-scene  system  to  support  ASW,  ASuW,  and  ISR  mission  planning  and  execution. 
Significant  collaborations  with  other  investigators  developing  and  integrating  UAS-based  data  and 
tools  will  be  also  leveraged.  Collaborations  with  Scripps  Institute  of  Oceanography  (SIO),  NSWC- 
Dahlgren,  and  other  PMW-120  EM  investigators  will  be  pursued  and  help  guide  our  sampling 
strategies,  assure  the  UAS  observations  taken  are  in  accord  with  NRL  goals,  and  enhance  synergistic 
development  of  UAS-based  TDAs.  Six  specific  technologies  that  will  be  developed,  tested,  and 
integrated  in  this  project  are  i)  tactical-scale  data  assimilation,  ii)  adjoint-based  observation  impact,  iii) 
air-ocean  coupled  modeling,  iv)  EM  tactical  decision  aid  systems,  v)  software  infrastructure,  and  vi) 
validation  and  verification  of  the  end-to-end  system. 

WORK  COMPLETED 

Work  that  has  been  completed  includes:  i)  development  of  UAS  Scan  Eagle  (SE)  data  processing 
software  (vl.O),  ii)  initial  testing/implementation  of  SE  data  processing,  iii)  evaluation  of  pre-TW13 
historical  COAMPS  forecasts,  iv)  successful  support  of  TW13  field  program  (including  3DVAR 
assimilation  of  SE  data  in  real  time,  development  and  evaluation  of  data  assimilation  innovation 
products  highlighting  differences  between  control  and  experiment  runs  using  SE  data,  initial 
evaluation  of  high-resolution  COAMPS  real-time  TW13  forecasts  and  adjoint  forecasts,  and  initial 
evaluation  of  SE  data  impact  on  Advanced  Refractive  Effects  Prediction  System  (AREPS)  output. 

i)  Scan  Eagle  data  processing  software  development :  Software  has  been  developed  that  includes  a 
new  SE  observation  decoder,  modification  of  original  UAS  time  averaging  routines  to  average  level 
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and  profile  observations  differently,  implementation  of  FGAT  (First  Guess  at  Appropriate  Time) 
software,  and  adjustment  of  data  thinning  parameters  (used  for  all  data). 

ii)  Scan  Eagle  data  processing  testing :  SE  UAS  data  processing  software  has  been  tested  for  an  initial 
sample  data  set,  and  extended  for  TW13  data.  Vertical  profiles  of  wind  speed  and  direction,  air 
temperature,  and  water  vapor  have  been  evaluated  for  the  lowest  ~1500-m  of  the  boundary  layer.  The 
sensitivity  of  vertical  M-profiles  to  both  slow-  and  fast-response  data  has  been  compared.  Results 
show  that  fast-response  data  captures  small-scale  structures  producing  critically  important  local 
trapping  regions  often  not  resolved  from  slow-response  data.  AREPS  has  been  evaluated  using  SE 
data  as  input  for  output  parameters  of  100  Mhz  VHF  Radio  frequency  and  X-band  radar. 

iii)  CO  AMPS  high-resolution  simulations :  Historical  COAMPS  simulations  for  July  2012  in  the 
TW13  field  program  area  (VAC APES)  have  been  conducted  and  analyzed  to  characterize  the 
climatological  environment  of  the  area.  Results  have  been  used  to  provide  guidance  for  model  setup 
in  preparation  for  the  field  phase  of  TW13. 

iv)  TW13  field  program  support :  Extensive  modeling  and  operational  support  was  provided  for  the 
July  2013  field  program  (see  Results  below). 

RESULTS 

a.  Scan  Eagle  data  processing  software  development 

A  significant  amount  of  software  development  was  made  to  the  COAMPS  data  assimilation  suite  to 
incorporate  UAS  data.  First,  UAS  observations  pass  through  aircraft  QC  procedures  originally 
developed  for  commercial  aircraft  observations.  The  UAS  observations  are  elevation  limited: 
observations  far  above  the  maximum  height  range  of  the  UAS  aircraft  are  automatically  rejected; 
observations  collected  at  15-m  elevation  or  lower  are  automatically  rejected  because  the  sensors  are 
not  shut  off  during  landing/waiting  for  take-off.  Second,  because  meteorological  sensors  and  GPS 
sensors  used  different  clocks,  synchronization  errors  can  produce  unusable  real-time  winds.  Thus, 
winds  are  not  used  (though  wind  errors  can  be  corrected  in  post-processing).  For  the  TW13  field 
program,  the  meteorological  sensor  clock  was  incorrectly  set  for  part  of  the  experiment  (off  by  hours), 
causing  all  meteorological  variables  to  be  rejected  because  they  were  too  “old”.  This  time  period  was 
-00  UTC  15  July  to  00  UTC  16  July.  Finally,  NAVDAS  performs  additional  data  thinning  (to  remain 
computationally  efficient),  thinning  data  in  areas  where  there  are  too  many  observations.  So  not  all  QC 
accepted  SE  data  are  used.  Figure  1  describes  the  new  data  assimilation  technique  called  First  Guess 
at  Appropriate  Time  (FGAT)  that  has  been  developed  by  NRL-MRY  and  tested  and  used  with  SE  data 
for  the  TW 1 3  field  experiment.  FGAT  is  a  technique  that  allows  data  to  be  used  at  the  most 
appropriate  time  within  the  data  assimilation  window  for  NAVDAS.  The  current  FGAT 
implementation  uses  piecewise  linear  interpolation  between  the  3-,  6-,  and  9-h  forecast  background 
fields.  Results  for  the  TW13  field  experiment  indicate  that  using  FGAT  provides  improved  analyses. 

b.  Scan  Eagle  data  processing  testing 

Figure  2  shows  the  impact  of  using  SE  data  in  COAMPS  analyses  for  the  TW13  field  program.  Two 
parallel  COAMPS  simulations  were  conducted  in  real-time,  one  assimilating  all  routinely  available 
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data  plus  SE  data  (called  twl3se)  and  one  assimilating  only  routine  data  (twl3).  Shown  is  the  1000- 
hPa  air  temperature  difference  between  two  COAMPS  nest-4  (1.33-km)  analyses  valid  12  UTC  16  July 
2013  for  TW13  (twl3se  minus  twl3).  The  triangles  are  locations  of  SE  observations,  with  the 
innovations  color  coded  the  same  as  the  temperature  differences.  Results  indicate  that  NAVDAS 
makes  the  analysis  using  SE  data  (twl3se)  warmer  in  the  vicinity  of  the  SE  observations  at  this  level 
(1000-hPa)  for  this  analysis  time.  This  is  in  agreement  with  the  positive  analysis  temperature 
difference  (twl3se  -  twl3).  Differences  further  away  from  SE  observation  locations  are  the  result  of 
accumulating  differences  between  analyses  (and  model  background  fields  based  on  the  6-h  forecasts) 
over  the  course  of  the  experiment. 

c.  COAMPS  adjoint  sensitivity 

An  adjoint  allows  one  to  compute  the  sensitivity  of  a  particular  forecast  output  (e.g.  the  modified 
refractivity  in  the  lowest  200  m)  to  changes  in  the  initial  state  in  a  mathematically  rigorous  and 
computationally  feasible  manner.  The  COAMPS  adjoint  identifies  regions  of  greatest  sensitivity  in  the 
initial  state  and  highlights  the  areas  where  additional  observations  (targeting)  could  improve  the 
forecast  the  most.  In  preparation  for  the  TW13  field  campaign,  the  first-ever  COAMPS  adjoint 
simulations  using  the  vertical  gradient  of  modified  refractivity  (M)  were  conducted  and  evaluated  for 
the  TW13  experiment  area.  We  have  computed  the  adjoint  sensitivity  response  function  of  the  6-h 
forecast  “M-profile”  modified  refractivity  to  initial  time  vertically  integrated  kinetic  energy  as  well  as 
initial  time  sea  surface  temperature  for  various  test  periods.  Results  from  this  example  shows  how 
forecasted  refractivity  in  an  area  of  interest  projects  back  to  where  SST  or  wind  speed  observations 
should  be  taken,  as  well  as  the  UAS  proposed  flight  profile  within  range  from  R/V  Knorr.  Using  this 
knowledge  we  employed  the  COAMPS  nested  adjoint  in  real-time  for  the  TW13  field  program.  Figure 
3a  shows  the  COAMPS  nest-4  (1.33-km)  adjoint  sensitivity  response  function  of  the  3-h  forecast 
vertical  gradient  of  modified  refractivity  (red  box  IVO  R/V  Knorr)  to  initial  time  vertically  integrated 
total  energy  for  00  UTC  16  July  2013.  Regions  of  large  sensitivity  occur  in  the  vicinity  of  the  Knorr, 
and  initial  condition  sensitivity  is  a  maximum  in  the  MABL.  The  variable  of  greatest  sensitivity  is 
moisture  (most  important  variable  to  observe).  If  the  optimal  perturbation  was  included,  there  could 
be  a  large  impact  on  EM  propagation,  as  illustrated  in  the  modified  refractivity  profiles  valid  09  UTC 
16  July  (Fig.  3b).  Note  that  the  profile  with  the  adjoint  optimal  interpolation  (red  line)  indicates  a 
ducting  regime  that  does  not  occur  for  the  control  (without  perturbation). 

d.  SE  data  impact  on  AREPS  output 

Currently  there  exist  several  EM  TDAs.  AREPS  is  the  Navy’s  sensor  and  communications  prediction 
software  system,  providing  radar  probability  of  detection  and  communications  assessments  for 
onboard  sensors  operating  in  coastal,  ocean  and  littoral  environments.  Refractive  effects  from  these 
mesoscale-rich  environments  are  introduced  through  ingest  of  NWP  model  forecast  fields  yielding 
vertical  modified  index  of  refraction  (M)-profiles  and  propagation  condition  summaries  along  with  the 
height-range  propagation  loss  for  multiple  bearings.  Figure  4  shows  an  example  from  the  TW13  field 
experiment  illustrating  the  impact  of  using  SE  data  on  EM  propagation.  Figure  4a  shows  a  map  of  the 
location  of  transmitter  (WYSJ-TV  Yorktown)  and  R/V  Knorr  used  for  calculations  of  EM  propagation. 
Point  propagation  summaries  for  the  Knorr  balloon  launch  valid  0933  UTC  16  July  and  the  two  18 
UTC  15  July  15-h  COAMPS  nest-4  (1.33-km)  forecasts  (twl3  and  twl3se)  are  shown  in  Fig.  4b. 
Results  indicate  similar  duct  thicknesses  between  the  observed  balloon  launch  and  the  COAMPS 
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forecasts.  Note  that  the  magenta  region  for  the  balloon  and  the  red  regions  for  the  COAMPS 
simulations  show  good  agreement  (thickness  ~200-m).  The  range  dependent  propagation  summaries 
are  shown  in  Fig.  4c  along  the  WYSJ-Knorr  path,  continuing  offshore  for  ~250-km.  The  forecast  for 
twl3  indicates  a  coastal  trapping  layer  that  increases  with  height  (indicated  by  red  line  in  top  panel) 
that  is  more  prominent  than  for  twl3se  (bottom  panel).  A  second  trapping  layer  exists  for  both 
forecasts  at  ~1 10-km  downrange  to  ~250-km.  The  associated  propagation  loss  coverage  diagrams  for 
the  three  vertical  profiles  along  the  WYSJ-Knorr  path  is  shown  in  Fig.  5.  There  is  similar  agreement 
between  the  balloon  and  both  COAMPS  forecasts  from  WYSJ  to  Knorr  downrange  to  ~60-km,  but 
more  loss  at  the  Knorr  in  the  lowest  ~100-m  for  twl3se  and  more  loss  beyond  Knorr  for  twl3  (>  ~100- 
km)  (as  indicated  by  the  dashed  white  ovals). 

RELATED  PROJECTS 

COAMPS  will  be  used  in  related  6.1  projects  within  PE  0601 153N  that  include  studies  of  air-ocean 
coupling  and  boundary  layer  studies,  and  in  related  6.2  projects  within  PE  060243 5N  that  focus  on  the 
development  of  the  atmospheric  components  (QC,  analysis,  initialization,  and  forecast  model)  of 
COAMPS. 


First  Guess  at  Appropriate  Time  (FGAT) 

•  Standard  3DVar  equation: 

*„  =  *(,  +  PbHT  [HPbHr  +  Pol  1  [y0  -  H(*„)] 

^  V  ^ 

Relative  to 

forecast  hour  3  6  9 

for  background  ^  ( 

•  Current  FGAT  implementation 
uses  piecewise  linear 
interpolation  between  the  3,  6,  9  h 
forecast  background  fields 

•  Runs  without  FGAT  compute 
innovations  using  6  h  forecast 
as  background 

•  FGAT  can  also  use  Lagrange 
interpolation  between  3,  6, 
and  9  h  forecasts 

Relative  to  g  q 

analysis  time 

Figure  1.  Description  of  the  new  First  Guess  at  Appropriate  Time  (FGAT)  technique. 
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1000-hPa  air  temperature  difference 


Figure  2 .  1000-hPa  air  temperature  difference  between  two  CO  AMPS  nest-4  (1.33-km)  analyses 
valid  12  UTC 16  July  2013  for  TW13  (COAMPS  with  SE  data  minus  COAMPS  without  SE  data). 
The  triangles  are  locations  of  SE  observations ,  with  the  innovations  color  coded  the  same  as  the 

temperature  differences. 


Figure  3.  a)  CO  AMPS  nest-4  (1.33-km)  adjoint  sensitivity  response  function  of  the  3-h  forecast 
vertical  gradient  of  modified  refr activity  (red  box  IVO  R/VKnorr)  to  initial  time  vertically 
integrated  total  energy  for  00  UTC  16  July  2013 ,  and  b)  modified  refr  activity  profiles  valid  09  UTC 
16  July  for  the  control  (black  line ,  without  perturbation)  and  with  the  adjoint  optimal  interpolation 

(red  line). 
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Figure  4 .  a)  Map  of  location  of  transmitter  (WYSJ-TV  Yorktown)  and  R/VKnorr  used  for 
calculations  of  EM  propagation,  b)  point  propagation  summaries  for  the  Knorr  balloon  launch  valid 
0933  UTC 16  July  (left)  and  the  two  18  UTC 15  July  15-h  COAMPS  nest-4  (1.33-km)  forecasts 
(twl3  and  twl3se),  and  c)  range  dependent  propagation  summaries  along  the  WYSJ-Knorr  path, 

continuing  offshore  for  ~250-km. 


Figure  5.  Propagation  loss  diagrams  for  the  Knorr  balloon  launch  (left)  as  described  in  Fig.  4,  and 

the  COAMPS  forecasts  (twl3  and  twl3se). 
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